function FF=mill_force_forward(UU)

n=UU(1)./2./pi.*60;
Vsup=UU(2);
b=UU(3);
t=UU(4);

%b=10;       %tool working lenght, mm
%Vsup=100;   %support speed , mm./s
phisup=0;   %support angle (???), grad
delta=60;   %cutting angle, grad
D=10;       %tool diam, mm
%t=1;        %mill depth, mm
%n=5000;     %tool rot. speed., 1./min
Z=3;        %blade count
T=120;      %total cutting time, min
alam=0.75;  %friction intensity koef.

Vs=Vsup./1000.*60;
R=D./2;
Sz=1000.*Vs./n./Z;
phiout=acosd((R-t)./R);
phicut=phiout./2;
acut=Sz.*sind(phicut);

v=pi.*D.*n./60./1000;
phiin=phisup+phicut;

p=1.9+0.04.*phiin;
k=(0.25+0.005.*phiin).*delta+(0.008+0.0018.*phiin).*v-(7+0.18.*phiin);

e=0.0008;
l=pi.*D.*phiout./360;
drho=e.*l.*n.*T./1000;
rho0=5;
arho=1+0.2.*drho./rho0;
if b>0
    if acut<0.1
        Fx1=(arho-0.8).*p+(8.*p+k+alam.*t./b).*acut;
    else
        Fx1=arho.*p+(k+alam.*t./b).*acut;
    end
else
    Fx1=0;
end


Fx=Fx1.*b;
Fxround=Fx.*l.*Z./pi./D;
Pp=Fxround.*v;
M=Fxround*R./1000;

FF=[M Fxround Fx Pp];

